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1 Pairs of quadratic hamiltonians 



In papers [H El [3l HI El [H [7] the problem of the commuting pairs of Hamiltonians quadratic 
in momenta was considered. 

Consider pair of Hamiltonians in the form 

H = apl + 2bpip2 + cpl + dpi + ep2 + /, (1.1) 

K = Apj + 2BpiP2 + Cpl + Dpi + Ep2 + F, (1.2) 

commuting with respect to standart poisson bracket {paiflis] = ^ap- The coefficients in 
formulas f ll.ip . (ll.2p - some (locally) analitical functions of the variables gi,g2- 

Theorem 1. Any pairs of commuting Hamiltonians (ll.ip - (ll.2p can he canonically transformed 

by 



osi osq2 
to the pair of the form 



H=^'-^\ K='-^^-'^^\ (1.3) 

Si — S2 Si — $2 



where 



TT C r ^ p2 , VSi{Si)S2{s2)Zs, Si{Si)Z 

(Sl - 32) 4(Si - SiV 



(1.4) 



U2 = S2is2)Pi - (,^_,^) ^1 - + ^.(^1, s.) 



Si - S2 
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V2 = 7v/:^9,, (y^S^i^i))^^) + f2{s2 
I \ S2 — Si J 



(1.5) 



for some functions Z{si, S2), Si{si) and fi{si). Poisson bracket {H, K} equals to zero if and 
only if 

and 

Proof We introduce new coordinates si,S2, such that the quadratic parts of H,K p.l|1.2p 
are diagonal: Let si, S2 be the roots of equations 

*(s, qi, q2) = {B- bsf -{A- as) (C - cs) = 0, (1.8) 
Then the canonical transformation 

(gi,g2,Pi,P2)-(si,s2,Pi,P2):pi = -(|fPi + |ip2), P2 = -(|f Pi + |ip2), (1.9) 



where = gi, ^2) under conditions {H, K} = transforms pairs p.ip .( fTT2ll to the form 

H=— K = ^-^ —, 1.10 

Sl - S2 Si - S2 

where 

Ui = Si{si)Pl + dPi + eP2 + f, U2 = S2{si)Pi + DP1 + EP2 + F, (1.11) 
where 

= J^2(i^'^ - + 2(fe. - + (^^« - CWJ') (1.12) 

We calculate a Poisson bracket between H and K. Then the coefficient of Pf, P2, P1P2 
equal to zero iff 



7 r dF{suS2) . ^/Sijsi) 32(32) Zs, ~ ^/Sl{sl)S2is2)Zs^ ~ dF{si,S2) 

where Z{si,S2), F{si,S2) - some functions, and 

■^^1,82 = Tzr- T (1-13) 

2{S2 - Sl) 

We apply the canonical transformation 

p^_p^ I dFisi,S2) ^ p^^p^ , dF{si,S2) 



OSi OS2 

to equate c?, E to zero. Then the coefficient of Pi, P2 equal to zero iff Ui, U2 have the form 
as in formulation of Theorem 1. And finally the free coefficient in Poisson bracket equals to 
zero iff the equation (I1.7P of the Theorem 1 is fulfilled just as expected. 

The general analytical solution of Euler - Darboux equation (II. 6p has near the line of 
singularities x = y the following expansion: 

00 00 
Z{x,y)=A + \n{x-y)B, A = ^ a,(x + y) (x - 2/)2\ B = Y^H^ + v) {x - yf\ 



where Oq and ai - some function. The other coefficients can be expressed by these two 
functions and their derivatives. For example, Bq = |aQ. 

We insert this expantion into p. 71) to obtain i? = 0. It is easy to check that any solution 
of the equation (11.61) with 5 = has the form 

Z{x, y)=zo + Six + y) + ix- yf ~ vf"^ (1-14) 

where g{x) - some function and Zq, S - some constants. We call the function g{x) as generating 
function for (I1.14p . Without the loss of generality we choose Zq = 0. The parameter 6, is very 
important for classification of hamiltonians from Theorem 1. 
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We find all the functions Z, corresponding the rational generating functions g. Choosing 
g[x) = x", we obtain the infinite set of polynomial solutions Z*^") for (II .6^ . In particular 

g{x) = l ^ Z(0)(x,y) = (x-y)2 

g{x)=x Z'-^\x,y) = {x + y){x-y)'^, 

g(x) = x^ ^ = i {{x - yf + 4(x + yf) (x - yf. 

All set can be obtained by using 'creating' operator 

2 2 9 I , s 

acting on Z^^\ The rational functions g{x) = (x— create one more class of exact solution 
of equation (11.61) . For example 



The solution corresponding the poles of order n > 2, can be obtained by differentiating the 
last formula by parameter /x. Because function Z is linear by g we obtained the solution Z 
with rational generating function g{x) = J2i Cj^^* + ^ijix — /ij)"-'. 

Hypothesis 1. For all Hamiltonians f ll.3p -( |T77l) generating function g is rational and has 
the form g{x) = where P m S - some polynomials with degP < 5, degS* < 6. 

In papers [SI [6] the following solution of the system p. 61) . p. 71) was considered: 

6 

Z{x,y) = x + y, Si{x) = S2ix) =^Cix\ 

i=0 

3 1 ^ . 

fi{x) = f2{x) = --CqX^ - -CgX^ + ^ kiX\ 

where Ci,ki - some constants. A very important fact is that Clebsch top and so(4)-Schottky- 
Manakov top [8l[9l[l0] are the particular cases of this model [6]. In paper [6] a full solution of 
Hamilton - Jacobi equation of this model was obtained in the form of some kind of separation 
of variables on a non-hyperelliptic curve of genus 4. 



2 Universal solution of Hamilton-Jacobi equation 

Let H and K have the form (ll.3l) - (|1.5p . Consider system H = ei, K = 62, where - some 
constants. Let pi = Fi{x,y), p2 = F2{x,y) - be its solution. We use short notation x n y 
corresponding gi n q2. Jacobi's lemma gives that if {H, K} = 0, then ^ = To find an 
action S{x,y, 61,62), it is enough to solve the following system 

f S = F„ As = f,. 

OX oy 
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We rewrite the system H = ei, K = 62 in the form 

pl + ap2 + b = 0, pl + Api + B = 0, (2.15) 



where 



Zx IS2iy) ^ Zy / 5*1(2;) 



A 



x-y \ Si{xy x-y\j S2{y) 

Zl , K - eix + 62 „ Z^ V2- eiy + 62 

-D = — T7 TTT + 



A{x-yr S,{x) ' 4(x-y)2 S2iy) " 
It easy to find that 

2by + Aa^ + 2aA^ = 0, 2Aay + aAy + 2B^ = 0. (2.16) 
Using ( ll.Gp and (11. 7p , it is easy to obtain the following identity 

Ab^-aBy + 2A^b-2ayB = 0. (2.17) 

Using a standard technique of Lagrange resolvents (see f.e. [H]), we rewrite system ( 12.15^ 
to a system 

uv = ^aA, (2.18) 
b B 

Au^ + i-u'^v - A—uv^ - av^ = 0, (2.19) 

that is equivalent to the qubic equation on v? . Let {uk,Vk), k = 1,2,3 be the solutions of 
(EH]), dim such that 

u\ + u\ + u\ = —b, vf + V2+v'^ = —B 

1 1 

U1U2U3 = — a^A, V1V2V3 = — A^a. 
8 8 

Then, formulas 

Pi = U1 + U2 + U3, p2 = vi +V2 + fs; 

Pl = U3-U1- U2, P2 = V3 - Vi - V2\ 

Pl=U2-Ui- Us, P2 = V2 - Vi - V3] 

Pl=Ui-U2- Us, P2 = Vi - V2 - V3 

define four solutions of (I2.15p . Consider the first of them. 

Lemma 1. For i = 1,2,3 following equations are fullfiled ^ = 

Prove. Differentiating equations (I2.18p and (I2.19p on x and y, we find Uy and Vx as the 
functions on u and v. Then expressing v through u, we obtain that Uy = Vx is equivalent to 
identities fl2T6D and fl2T7ll . ■ 

Lemma 1 means, that in variables ui,U2,U3 we find "particular" separation variables. 
Really S = Si + S2 + S3, where S is the action, and functions Sj defined from a system 

d_ - 
OX ay 
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Let's 



2x — yyx — ^ 2x — yyy — ^ 

It easy to see that pair {u, v) for all ^ are a solution of (|2.18l) . If Z is a solution of (11.61) . then 

du dv_ 

dy dx ' 



du _ dv Using this fact we introduce a function a{x, y, ^) so that 



da da 

dx ' dy 

In a case of rational function g, corresponding function Z is expressed through quadratic 

radicals and the function a can be obtained. Let's Y = -—. 

d^ 

After multiplication of expression (I2.19P by expression 

ZxZy 

left side of (12.19p can be written in the form 

x-y\ 4:{x - Q{x - y)/ x-y\ A[y - e.)[x - y) / 

Proposition 1. Let the expression lll.6\) . ( f i. 71) be fulfilled . Then the expression ( \2.20h 
is a function of Y and C, variables only. 

Prove. We assign the function (12.201) as \E'(x, y,^). Consider Jacobian 

j_d'^dY d^ dY 
dx dy dy dx 

We change ^ and ^ to || and ||, respectevily, then Jacobian J equals to zero identically 
taking into account HIM . (II. 7p . ■ 

Due to Proposition 1, the relation '^{x,y,^) = can be rewritten in the form Y) = 0. 
One can find the function by assuming y = x. 

Equation Y) = defines a curve, and the differentials of this curve define the function 
of action S. 

We note ^k{x,y), where k = 1,2,3, the roots of cubic equation ?/,^) = 0. 
Theorem 2. The function of action S has the form 

3 

S{x,y) = Y,{ ^(^^ 2/, ^k) - / YiO ) , (2.21) 

k=i 

where Y{^) - alebraic function on the curve 0(^, Y) = 0. 
Prove. We obtain 

^33 3 

— S(a;, y) = ^ cT^ix, y, ^k) + ^{ (^d^, y, 6) - Y{^k)}^k,x = ^Uk=Pi- 

k=l k=l k=l 



Analogously 

—S{x,y)=p2. ■ 

3 Case of cubics 

Consider a case when the curve (12 .201) can be written in the form 4>{C,,ri) = 4^ 

so that points {C,i,i]i),{C,2,f]2),{C,3,f]3) lis on a straight line, that equivalent to definition 

rj = ^a{x, y) + h{x, y), its substitution into (p, gives a curve \E'(a;, y, C,) = 0. 

Formula (12.201) gives a curve in a new variables ^,?7 —62 + ei^ + = 0, where 

^) at a; ^ or y ^ 0. 

Using reversible curve equation Ci{^,ri) = ^ = /(^) using rj, we find the expressions 

for a{x,y), b{x,y) 

aix,y) = M^IM^ bix,y) = ^^^^^ ' "^^^^ 
X ~ y X ~ y 

On the other hand the equivalence of the curve F) = h 77) = gives 

YxVy = YyVx ^ U^Vy = ^^Tj^ ^ H - y)Z^r]y = (x - O^yVx, 

or Z = Z{a), = -ya^, by = -xay. 

4 Examples 

In this Section we consider all the pairs of Hamiltonians known at the moment f ll.3l) - p.7p . 

4.1 Class 1 

For the models of this class 

S^ = S2 = 5, /i = /2 = /. (4.22) 

Theorem 3. Let 

9 = t G = G-^S', f = -^-^G'-^-S", 

^ S 10 S3 12 ' 

where 

S{x) = s^x^ + S4X^ + SsX^ + S2X^ + S\X + Sq, G{x) = g^x^ + g2x'^ + gix + g^, 

where Si, gi, 6 - some constants. Then functions S, f and function Z, corresponding (see. §1) 
generation function g, fulfill the systems (11.61) . (11.71) . 
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Remark. Parameter S from Theorem 3 coinsides with parameter S from f ll.l4p . Consider 
the case 5 = in the formula (fTTiD . Then all pairs of Hamiltonians ^^-^1^, (14221) . that 
fullfil this condition are described by Theorem 3. 

Consider a general case 

S{x) = S5(x - - fl2){x - fi3){x - fi4)ix - yUs), 

where S5 7^ and all roots Hi of polynomial S are distinct, then the function Z has the form 

5 

y) = Yl -y), (4.23) 

i=l 

where z/j - some constants. Coeffitients gi and S are expressed through constants Uj from 
(I1.14p . For example, 26 = —J2 ^i- Function / is defined by 



fix) = -^Y.^^^ + hx + ko, 
lb ^-^ X — Ui 

1=1 



where ki,ko - some constants. 

Calculation for a function (14.230 gives 



(^{x, y, = - o 2^ log , , ^ , (4.24) 



4tr \i-^^^)^{^-ay-0' 

Algebraic curve has the form of hyperelliptic curve of genus = 2 

0(y, = 5(0>^' + /(O - + 62 = 

Steklov top on so(4) [12j is a particular case of Theorem 3. 
4.2 Class 2 

Functions Z for the models of this class are the special cases of the functions Z of Class 1. 
But this Class contains much more parameters them Theorem 3. 

Such functions Z can be defined as the solutions of system 

^x. = ^^ = ^f^(^)^.^., (4.25) 
where U - some functions of one variable. 



Remark. It easy to see that this class of solutions of Euler - Darboux equation Z 



— — - coincide with the class of solutions of the form 



xy 



2{y-x) 

Z = F 



hix)-hiy) \ 
x-y J 
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where F and h - some functions of one variable and U = F" / F'"^. 
Lemma. The system (14.25^ is compatible if and only if 

U=l^, B{Z) = b2Z^ + biZ + bo, 
/ Jd 

where bi - some constants. 
In a case deg B = 2 

Z{x, y) = {x - iii){y - Hi) + - /i2)(y - /i2), (4.26) 

where 62 = 1, &i = 0, 60 = -(/^i - /^2)^- 
If deg 5 = 1, then 

Z{x,y) = ^+-{x + y), (4.27) 



2 

61 = 1, 62 = &o = 0. 
If deg 5 = 0, then 

Z{x,y)=x + y. (4.28) 
1. Consider function Z of the form 04.260 . Then 

S{x) = {x- /ii)(x - /i2)P(x) + (x - Aii)^/2(x - /i2)^/^Q(x), degP < 3, degQ < 2, 
fix) = fo + /ix + fc2(x - /.0^/^(x - /.2)^/^ + ~ ^ "^^^'^ "^^^^ 



16 I X — /ii X — H2 

+ (X - /.0^/^(x - /.2)^/^ - 

32 I X — /ii X — yU.2 J 

In a case when Q = 0, A;2 = 0, These formulas coinside with corresponding formulas of Class 
1. The functions a, Y are defined the same formula (I4.24p as for Class 1 : 

2 2 / 

. iv^i - ^v?/ - + - Wx - /^i ^ 1 V(a;-Mi)(l/-^i) 

2^ Vx^VtM^ 4^ (e-/ii)V(x-0(z/- 

Algebraic curve in this case has the form 

[^«(0 + r/^/(e)]>^' - [kniO + rikiiO] = 0, (4.29) 



where 



Sr{x) = (x - /ii)(x - fi2)P{x), Si{x) = (X - /ii)(x - fl2)Q{x), 

16 I X — yUl X — /i2 J 

A:.(a:) = A:. - ^^.(Mi - ^.)^ - ^(Mi _ ^,)|^ _ ^ 
32 16 y. X — jii X — jjL2 



^ v/e^v/e^V 16(e-^i)2(e-/^2)2>^^' 



Expressing F as a function of //) and substituting to (14.290 . we obtain 10-parameter 
cubic in (^,?7), variables. 

So in a general case the curve (f>(Y,^) = 0, is a covering over an elliptic curve We obtain 

^ - Ail , ^ - At2 



V = -7E== ^== + 



yx-fii Vy-i-ii y'x-/i2 I Vy-t^2 ' 
\Jx-iii \Jx-ii\ \Jy-ii.\ 

therefore points (^2,'72), (^3,^73) lie on a straight line. 

2. For the function Z of the form (14.27^ we have 

S[x) = xP{x) + x^''^Q{x), degP < 3, degQ < 2, 

Ibx 61^ X 

^ + \fxJy 

The function Y is defined by F = ^ . The curve in this case can be written 

4ev/^v/^ 

in the form f l4.29p . where 

SrIx) = xP{x), Si{x) = xQ{x), 

1 1 

knix) = -62 + eix - /o - fiX + -— P(x), ki{x) = -r^Qix) - fg, 

Ibx Ibx 

V 



In rj) variables it also has the form of arbitrary cubic. Formula rj = gives the fact 

that points (^1, r/i), (^2, {(,3, Vs) lie on a straight line 
3. For the function Z, given by (14.281) . we obtain 

S{x) = S^X^ + S^X^ + S^x"^ + S3X^ + S2X^ + SiX + Sq, 

f{x) = -j^S"{x) - -^g(x) + hx' + hx + u 
4U SZy/x 

In this case Y = — ^ . Algebraic curve 

SiOY' - FiOY' - + ^S^'^iO - y) - II = 0' F(0 = -e2 + e,^ - /(O 

and in (^,77), variables where rj = ^'^ — -, has the form of arbitrary cubic. Because 

4y 2 

V = ^{x + y)- xy, points (6, ^71), (6, ^72), (6, V3) belong to a straight line. 
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4.3 Class 3 



We introduce 'non-symmetrical Hamiltonian (11.3p - p.7p if Si{x) 7^ 5*2 (x), or 7^ 

f2{x). 

Theorem 4. [7j In non-symmetrical case the functions Z, Si, fi is the solutions of I \1M . 
(11.71) if and only if 

5 = 0, g = ^, S,,2 = WH±MH'/\ f,^2 = -^T2MH-^/^±aH^'\ 

where g - generation function of Z, 

W{x) = WsX"^ + W2X^ + WiX + Wq, H{x) = h2X^ + hix + Kq, 

M(x) = m2X^ + rriix + rriQ. 
Here Wi, hi,mi,a - some constants. 

Consider the general case H{x) = {x — ^i){x — ^2)- Algebraic curve in this case is defined 

by 

^(e, Y) = -62 + - S^^^f S+ 

(4.30) 

where 

Substituting 

1 (/il-;U2)2?7 



Y 



4 (e - /i2)(e - /Xl)V^^(/i2-/il)-8(e-/Xi)(e-/X2) 

into (14.30p . We obtain the cubic in variables (^,77) with a full set of ten independent 
parameters. It easy to see that 77 = a{x, y)^ + h{x, y), where a,h - some functions. 

Therefore in the cases of Class 2 and 3 the algebraic curve is non-hyperelliptic covering 
over the elliptic curve. The dynamics of the points (^1, Yi), (^2, Y2), (^3, 13) on this curve (see. 
theorem 2) defines the following condition: their projection on the elliptic base (^1, r/i), (^2, 112)1 ('Cs; ^3) 
lies on the straight line. 

Hypothesis 2. Any pair of the Hamiltonians (ll.3p - (11.7p belongs to one of above three 
classes. 



5 Appendix 1. Steklov top 

We show that the case of Steklov top on so(4) is a particular case of Class 1 after restriction 
on the symplectic leafs. Hamiltonian and the additional integral in this case have the form 

= (Si, ASi) + (Si, BS2), K = (Si, ASi) + (Si, 5S2), 
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where 



A = —Q^ diag(—^, -^), B = a diag(ai, 02, as) 
ai a2 ai 



A = —diag{a\, al, al), B = a diag{ — , — , — ), 

Oil 012 «3 

a = aia20is. Here Sj - three-dimensional vectors with components S". It is easy to see that, 
H and K commute under a spin Poisson bracket 

It is convenient to chose n = —2i. 

We fix the Casimirs for spin bracket: (S^, S^) = j|, k = 1,2. Then the formulas 

Sk = TTkKiQk) + yK'(Qfc), W K(Q) = ((Q2 _ 1), ,(q2 ^ 2Q), 

define the Darboux coordinate tti,tt2, Qi, Q2 for the simplectic leaf of Poisson manifold 
with coordinatesSfc, k = 1,2. As this transformation is linear by momenta TTfc, as a result, 
we obtain a pair of commuting Hamiltonians quadratic in momenta under the bracket 
{T^ajQis} = Sai3- Cousidcr the canonical transformation that transforms their pair to the 
form (O-drS!), gM)- 

We apply the canonical transformation 

Pi = vriv/KQO, P2 = n2./R(^), dX = ^^, dY = -^^, 

Vr{Qi) V^iQv 

where 

riQ^) = (K(gi), AK{Q,)), R{Q2) = (K(g2), ^K(Q2)), 

to obtain 

jl(9'i{X) 3^g[{X). 



g,{X) 2'g,{X)' 



H = Pl + 2P1P2V + J2P1VY + J1P2VX + -Jlj2Vx,Y + J ( 
K = P^ + 2P,P2W + j2P,Wy + nP2Wx + IjiJ2Wx,Y + | - ^(^; 

2 6 \g2{Y) 2 g2{Y) 

where 

^ m.iB^m) ^^^^ ^ m^B^m) 



We apply the canonical transformation [Pi, P2, X,Y) {pi,p2,x,y) of the form 
dX=U^ + ^i^), dY = -U^^ + ^\. 



Pi 



ji + 32 fy\ fE^r^ ( , Ji + 32 fx' 

Pi 



12 



pi 



31+32 y 



4(x — y)\ X 



P2 + 



ji + 32 [x 



A{x-y)]l y 



to obtain (Ol-dllSl), KTh . where 

S{x) = — 4x(l + + alx){l + alx), Z{x, y) 



-Ji{x + y) - 32^/xy 



fix) 
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